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We present a relativistic chiral mean field (RCMF) model, which is a method for the 
proper treatment of pion-exchange interaction in the nuclear many-body problem. There the 
dominant term of the pionic correlation is expressed in two-particle two-hole (2p-2h) states 
with particle-holes having pionic quantum number, J". The charge-and-parity-projected 
, relativistic mean field (CPPRMF) model developed so far treats surface properties of pionic 

■ correlation in 2p-2h states with = 0~ (spherical ansatz). We extend the CPPRMF model 

by taking 2p-2h states with higher spin quantum numbers, J n = l + ,2~,3 + , ... to describe 
the full strength of the pionic correlation in the intermediate range (r > 0.5 fm). We apply 
the RCMF model to the 4 He nucleus as a pilot calculation for the study of medium and heavy 
nuclei. We study the behavior of energy convergence with the pionic quantum number, J" , 
qq 1 and find convergence around J^ax — 6~ . We include further the effect of the short-range 

repulsion in terms of the unitary correlation operator method (UCOM) for the central part 
of the pion-exchange interaction. The energy contribution of about 50% of the net two-body 
interaction comes from the tensor part and 20% comes from the spin-spin central part of the 
pion-exchange interaction. 

O 

§1. Introduction 

Chiral symmetry is the most important symmetry as a key to understanding the 
mechanism of mass generation and pion dynamics in the subatomic world governed 
by the strong interaction. 1 ' At the hadron level, chiral symmetry is described using 
[~». \ the linear a model introduced by Gell-Mann and Levy!^ The Nambu-Jona-Lasinio 

(NJL) model is used for the mass generation and chiral properties of quark matter.® 
Pion emerges with a zero mass field as a Nambu-Goldstone boson by spontaneous 



chiral symmetry breaking. With a small explicit symmetry breaking term, the pion 
acquires a small finite mass as compared with other mesons. Pion plays an important 
role as a mediator of intermediate- and long-range parts of nucleon-nucleon interac- 
tionPJ 1 It is extremely interesting and challenging to construct a nuclear system with 
the chiral Lagrangian, where the pion plays an important role in the nuclear struc- 
ture. We therefore study the property of a nuclear system with the linear a model 
of Gell-Mann and LevyP The proper treatment of the pion-exchange interaction 
in a way consistent with the many-body description becomes a key issue. For this 
purpose, it is necessary to construct a relativistic framework describing naturally the 
mass generation of hadrons and at the same time treat the pion-exchange interac- 
tion. We would like to develop a theoretical framework to handle the pion-exchange 
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interaction together with the partial restoration of chiral symmetry in the nuclear 
system. 

We have started this program of treating the chiral sigma model Lagrangian 
in terms of the relativistic mean field approximation by introducing parity-mixed 
intrinsic single-particle states.^'^ We further perform the projection of the parity 
and, in addition, the charge of the total wave function from the parity- and charge- 
symmetry-broken intrinsic nuclear wave function, which is named charge-and-parity- 
projected relativistic mean field (CPPRMF) modelP® We have shown that the 
pion plays an important role in generating an effect of spin-orbit splitting and leads 
to jj-closed shell magic numbers.^® This charge and parity projection was also 
performed in the nonrelativistic Hartree-Fock schemeP®^^ 

The CPPRMF model provides a specific feature that the 2-particle states of 2p- 
2h component have a spatially compact structure due to the pseudo-scalar nature 
of the pion exchange interaction. The Op-Oh to 2p-2h matrix elements of the pion- 
exchange interaction to provide this compact structure, and the energy variation with 
respect to the spatial wave function in the CPPRMF is extremely important. This is 
a critical difference of the CPPRMF model from any perturbative treatment of pion- 
exchange interaction.^'^'™ In the CPPRMF model, the dominant contribution 
of the pion is expressed in terms of 2p-2h configurations with particle-hole spin- 
parity to be J n = CP (spherical ansatz). The pionic energy contribution per nucleon, 
however, decreases rapidly with mass number in the systematic calculation of various 
medium and heavy nuclei. ®' This result seems inconsistent with that of the chiral 
perturbation theory for nuclear matter by Kaiser et alP^ They show that an iterated 
one pion-exchange of Hartree term gives a large energy contribution to the binding 
energy of nuclear matter. Our energy systematics of the CPPRMF model also seems 
to be not accordance with a result of the variational Monte-Carlo (VMC) calculation 
of light nuclei obtained by the Argonne-Illinois group This fact shows that the 
pion-exchange interaction has a volume effect. Hence, it is an important next step 
to include the volume effect in the CPPRMF model. 

In nuclear matter, the momentum is the good quantum number and, particularly, 
the pseudo-scalar nature of the pion exchange interaction produces the contribution 
of large momentum component. If we were to take into account the volume effect 
of the pion exchange interaction, we have to take angular momenta of the pion, L-p, 
up to the value satisfying q ma , x R = Lp ax , where R represents the nuclear radius and 
<7max represents the maximum momentum of an exchanged pion in pion-exchange 
interaction. If we take the maximum momentum as twice of kp ~1.4 fm -1 p3^D3 
we ought to take the angular momentum of the pion up to L^ ax ~ 6 for 4 He. 
Hence, we have to extend the CPPRMF model to include higher angular momentum 
components of the pion. The angular momentum of the pion, Lp, and the total 
spin-parity quantum number of the particle-hole states, J w , are related as CP ® 

Up ' = J 71 ", where CP is an intrinsic spin-parity of pions. Owing to this relation 
between both quantities, we call the total spin-parity of the particle-hole states, J 71 ", 
"pionic quantum number" hereafter. The pseudo-scalar nature of the pion-exchange 
interaction is expressed using the 2p-2h components with pionic quantum number, 
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J w = , in the CPPRMF model so far. As stated above, it is natural to extend the 
CPPRMF model to include 2p-2h components with higher pionic quantum number; 



In this paper, we formulate the relativistic chiral mean field (RCMF) model to 
treat full strength the pionic correlations in the nuclear many-body system. The 
RCMF model is based on the mean field picture, where all the nucleons move in 
the mean field generated by other nucleons. The pion has the pseudo-scalar nature, 
which is then treated in 2p-2h components with pionic quantum number; J 71 " = 
CP, 1 + , 2 _ , 3 + , .... Although the probability of nucleon pairs being brought up to 
high momentum states is of the order of 10 ~ 15%, the energy gain associated with 
the pion exchange interaction is essential to bind nucleons in the nucleus. 

This article is arranged as follows. In §2, we discuss the formulation of the RCMF 
model for finite nuclei. In §3, we apply the RCMF model to the 4 He nucleus as a 
pilot calculation. The numerical results are presented in this section and detailed 
discussions are made for the properties of the resulting wave functions and energy 
components. A summary and the outlook are presented in §4. 



We develop here the RCMF model for the description of finite nuclei. The RCMF 
model is a natural extension of the CPPRMF model, which properly takes into 
account the pseudo-scalar nature of the pion-exchange interaction. As discussed in 
the introduction, we include all possible 2p-2h configurations to take into account the 
full strength of the pionic correlations. We take the pion-exchange interaction with 
the form factor due to the finite size effect of the nucleon. We calculate explicitly the 
direct part of the pion-exchange interaction and drop exchange terms as the natural 
extension of the CPPRMF model.®'® 

2.1. Chiral model Lagrangian 

We start with the linear sigma model Lagrangian of Gell-Mann and Levy, where 
the pion field appears symmetrically with the a field. The pseudo-scalar pion-nucleon 
coupling in the linear a model leads to an unrealistically large attractive contribution 
through the strong coupling between positive and negative energy states, and we 
have to treat properly the effect of the negative energy states. We thus employ 
the nonlinear realization of chiral Lagrangian, which is obtained by the Weinberg 
transformation of the linear a modelP® We take the lowest order term in the pion 
field, and the Lagrangian density in the nonlinear representation is written as^®'® 



r = 



0-,l+,2- 3+ .... 



§2. Formulation 




(2-1) 



where 



= Tp{i^ -M-g a a- g^co^ 
+ \d^a - ^m a 2 a 2 - A/^cr 3 - '- 
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and 

C n = -^^aWty + ~d^ a dV - ^m n 2 7T a2 . (2-3) 

The nucleon and ui meson masses are obtained by a meson condensation in the 
vacuum. The effective nucleon mass and the effective uj meson mass are given by 
M* = M + g a a, and m* = m w + g^a. Masses and coupling constants are set as 
M = g a fw, "i 2 = A* 2 + A/ 2 , m 2 = n 2 + 3A/ 2 , and m u = gZf n . We take the empirical 
values for masses and the pion decay constant as M = 939 MeV, m w = 783 MeV, 
= 139 MeV, and f n = 93 MeV. The cr-nucleon coupling constant, g a , and 
the a-uj coupling constant, g^, are fixed by the relations, g a = M/f n = 10.1 and 
9~uj = m uj/fir = 8.42. The strengths of the a meson self-energy terms depend on the 
a meson mass, m a , through the relation A = (m 2 — m 2 )/2f 2 . The a meson mass 
and w-nucleon coupling constant, g w , are the free parameters. We introduce the 
pion-nucleon axial vector coupling constant, g^, in this Lagrangian. We set gA = 
1.25, which is related with the coupling strength in the free-space uNN scattering 
by the Goldberger-Treiman relation.^ 

2.2. Ground state wave function 

The total wave function of the whole nuclear system is written as 

= ® ■ (2-4) 

The meson wave function, \&m, is a coherent state for the a and oj mesons. As for the 
nucleon wave function, we write the projected wave function of CPPRMF model, 

& = N \0p-0h) + Y,N l \2 P -2h) l + Y,N k \4p-4h) k + --- . (2-5) 

i k 

The contribution of the pionic correlation is dominantly expressed in terms of 2p-2h 
states, where particle-holes of 2p-2h states have J n = CP. For 4 He, the probability 
of 4p-4h states is negligibly small.^ 1 Hence, we neglect the configurations with more 
than 4p-4h states. We introduce higher multipoles by coupling particle-holes of 2p-2h 
states to various spins and parities with pion quantum number. It is a straightforward 
extension to introduce higher mutipoles for the pionic quantum number, J 71 ". Hence, 
we would like to write the total wave function for the nuclear part, <Fat, as Op-Oh 
state and the combination of 2p-2h states, 

\V N ) = a \0p - 0h;b ) + ^ai\2 Pi - 2h i; bi), (2-6) 

i 

where i denotes various 2p-2h states. The normalization condition of the total wave 
function is (^nI^n) = ao^o + J2i a i a i = 1- The amplitude, «j, represents the 
weights of 2p-2h states. They are variational parameters for total energy minimiza- 
tion. The Op-Oh wave function, \0p — Oh;bo), represents the relativistic mean field 
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Fig. 1. 2p-2h Hartree diagram. 

(RMF) ground state, which is obtained by solving the RMF equations for Hamilto- 
nian, H a ^. Nucleons occupy mean field levels up to the Fermi surface. The 2p-2h 
wave function is constructed as 

r / / / i (oo) 

\2pi - 2hi- hi) = \p-h; JMTM T {bi)) ® \p' - h';J M T M^h)) . (2-7) 

L J i 

A single-particle wave function for 2p-2h states is written as 

where G nft (r) and F nft (r) represent the radial parts of the upper and lower compo- 
nents of the Dirac spinor. The orbital angular momentum and spin state are given as 
yKm(r,a) = Y^ mi ,m 3 ( lm i 1 / 2m s\j m ) Y im l (r)Xrn s (o-). The label n distinguishes radial 
wave functions of states with the same k, m quantum number. The label k defines I 
and j as k = — (Z + l) for I = j — 1/2 and k = I for / = j + 1/2, and R = —k represents 
the opposite parity state. C( r ) denotes isospin state. We follow conventions as in 
Refill]). 

The upper radial part is given by Gaussian function as, 



G nK =MN(l,b t )r l e W (-^), N(IM) = J ^sLj+sy ^ 



2 



where 6« denotes a range of Gaussian function and M is the normalization factor of 
the spinor, ip n j mT . In the case that the particle state has the same quantum number 
as the hole state, we take the next form to make the particle state orthogonal to the 
hole state for the upper radial part, 



2 

G nK = (M +N 2 r 2 )N{lM)r l ^v(-^ I ), (2-10) 

where M\ and A2 are fixed by the normalization and orthogonalization with the hole 
state. As for the lower radial part, we take the following form, 

1 / d k + 1\ ^ , 
F ™ = 7^7{-r + ) G nK , (2-11) 



2M\dr r 

which is the relation of the lower radial component with the upper radial component 
for plain wave function for small momentum, E ~ M. We take this assumption to 
minimize the number of variational parameters. 
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2.3. Hamiltonian 

In the relativistic mean field approximation, the source term of the pion disap- 
pears for parity-conserved single-particle states, and therefore, the solution obtained 
at the mean field level does not include the effect of the pion. We take this relativis- 
tic mean field ground state as an initial state for the RMF Op-Oh ground state. We 
express 2p-2h states in terms of Gaussian functions with various range parameters 
for pionic correlations. 

The relation between the Hamiltonian density and the Lagrangian density is 
written as 

H = y^<P-C, (2-12) 

where <j> denotes the nucleon field, tp, and ir, a, oj meson fields. The Hamiltonian 
density for a and oj fields is 

n a ,u = • V + M + g a a + g^uj^ (2-13) 

+ Iva • W + \mla 2 + Xf w a 3 + ^a 4 

- ^Voo° • Voo° + -(V x Co) • (V x Co) - -mlu^ 

- giU^^ - hjla 2 ^^, 

where we consider the static case, and the time derivative terms, docfr, vanish. We 
take the mean field approximation for a and oj field operators, 

a — >(&\<r\V)=a, (2-14) 
co — > {V\oj^) = S^ouo. (2-15) 

The expectation value of the spatial component, {Co), vanishes under the rotational 
symmetry condition. Hereafter, we write just V and 'w' as the classical fields, a 
and Oo, respectively. The Hamiltonian for a and oj mesons is obtained as 

H a , u = J d 3 xH^ (2-16) 

= J d 3 x^(-ia- V + jo(M + g a a) + g u oj)ip + £ meson (a,oj)} , 

where 

£mcson(c7,0j) = ^(-V 2 + ir&G 2 + A/^ 3 + (2-17) 
1/ ^7 2 1 2\ 2 ~2 t 2 1—2 22 

--(-V +mJoj -g u f w au - -g w a oj . 
The Hamiltonian density for the pion part is given by 

= 3- ^757 • Vvr a r> + \vit a • Vvr a + \m n 2 itV. (2-18) 
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Using the Eular-Lagrange equation, 

„ dC dC 



0, (2-19) 



we obtain the Klein-Gordon equation for the pion field as 

(-V 2 + V = |f V • ^757r a V, (2-20) 

where a = 0, ±. From this Klein-Gordon equation, we obtain the integral form for 
the pion field, 

— Tfi-K | X' — y\ 

d3 y ^—f\ V y ■ m l5 lT a ^{y) (2-21) 

By using this integral form and the Klein-Gordon equation of the pion field, the 
Hamiltonian for the pion-exchange interaction is written as 

H n = J d z xU n (2-22) 
1 f d 3 q 



2 J (2ttY 



d 3 xd 3 yip(x)-f5 ; y ■ qT a ip(x) 
2/,r cp + m n 



9 A \ 2 e^ s ~^> l - 



The total Hamiltonian is given by 

£ = + H v . (2-23) 

2.4. Pionic energy minimization 

We take <F as the trial wave function with variational parameters, c^, 6«, cr, and 
uj. As for the nucleon part, we search for the energy minimum within the (ojj, bi- 
plane. For this purpose, we solve the following two-energy minimization equations 
simultaneously One is given by 

d 

— mH-Em = 0, (2-24) 

OCii 

and the other is given by 

±(q,\H-E\V} = 0. (2-25) 

The first energy minimization condition provides the variational parameter «j. The 
second one is very important from the viewpoint of the pseudo-scalar nature of the 
pion-exchange interaction. The finding obtained by CPPRMF method is that the 
particle states of 2p-2h states have compact distributions due to the high-momentum 
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components caused by the pseudo-scalar nature. In the second step of minimization, 
we can achieve the full strength of pionic correlations by using a small number of 
bases. We make "pionic optimization" in this size variation. 

We obtain equations for the meson mean fields by varying the total energy with 
respect to a and ui classical fields as 

-^(V\H\g,) = 0, (2-26) 

where eft means a and uj. Differential equations for the a and co meson fields are 

(-V 2 + ml)a = -g aPs - 3XUa 2 - Xa 3 + gZ 2 f^ 2 + g^au 2 , (2-27) 

— * 2 2 2 2 2 

(-V + m u )U) = g u p v - 2g 0J f^ato - g u a u, (2-28) 

where the scalar and vector densities are calculated using the Op-Oh and 2p-2h wave 
functions presented in Eq. (2.6) as 

Ps = {¥ N \: W : |<%>, (2-29) 

Pv = {& N \-- V'ToV' : \&n)- 

We solve these four equations given by energy functional variations, Eqs. (2.24), 
(2.25), (2.27), (2.28), (2.29), self-consistently. We note here that these meson fields 
are affected by the 2p-2h pionic correction terms through each density shown in Eq. 
(2.29). The Op-Oh RMF ground state is re-decided by these a and uj meson fields for 
every iterative calculation until we obtain the self-consistent solution. 

2.5. Matrix element of the pion- exchange interaction 

We show how to calculate the matrix element, (2pi — 2hi;bi\H n \0p — Oh;bo), 
in this subsection. At first, we consider particle-hole states induced by the pionic 
response with quantum number, JMTMt, and two particle-hole states couple to + 
ground state. We put the suffix, "p" , for the quantum numbers of the particle state 
like 

n p>jp> Ipi an d we put the suffix, "/i", for the quantum numbers of the hole state 
like rih,jh^h- lp an d lh represent the single-particle orbital angular momentum of 
the upper component of particle and hole states, respectively. The lower components 
are represented by L and 1^, respectively. The pion carries momentum q between 
two nucleons. Following the method shown in Ref. I22|) . we write down the matrix 
elements of the particle-hole excitation induced by the pionic response as follows, 

(p - h; JMTM T |7o757 • qr x exp(-iq ■ f)\6) (2-30) 
= iq5 T \ Sm t X { ((lp 2 )jp ( l h g )jh 5 JM \a ■ gexp (-iq ■ r) 1 6) 
+ {( l p\)lp( l h]j)3h-, JM\a ■ qexp(-iq- r)|6)| 

= iq5 T ^ MT ^Y*j M {q){Y, a LJ F it{<l) +Y, a JLFi L *(q)}> 

L L 

where for the upper component, 

F^(q) = 2(LS; J\j p j h ; J)(-i) L ^4n(2l h + l)(Z ft 0L0|Z p 0) (2-31) 
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x / r drj L (qr)G npjplp (r)G nhjh i h (r), 



(LS; J\j p j h ; J) = [(2L + 1)(25 + l)(2j p + l)(2j h + 1) 

l p lh L 
1 1 S 



1/2 



2 2 



a JL = (-1)(J010|L0) = { 



2J+1 



1/2 



V2J+1 J 



1/2 



for L = J - 1, 
for L = J + 1. 



and for the lower component, 

F^(q) = 2(LS; J\j p j h ; J)(-i) L J 4*{2l h + 1)(40; Z0|f p 0) 



(2-32) 



X ^ r dr j z {qr)F npjpl - p (r)F n ^- (r) , 



(LS; J|jpj fc ; J) = (2L + 1)(25 + l)(2j p + l)(2j A + 1) 

> 



nl/2 



X < 



Z p t/j L 
lie 

„ jp jh J > 



a jZ = (-1)(J010|L0) 



2J+1 



^2J+l ) 



1/2 
1/2 



for L = J - 1, 
for Z = J + 1. 



|0) represents the Op-Oh RMF ground state. Using this pionlike response matrix 
element, we construct the matrix element of the direct term for the correlation 
between Op-Oh and 2p-2h states as follows, 

1 r d 3 q 



r _ < 9A A2 ex P(^' (n ~ ^2)) 



m 2 + g 2 



(2p - 2/1(70757 • (Fa 
= K" ( 2T )2 ) E( JMJ/M 'I° ) E (™ T T'il4|00)£(-l 



70757 • QT\\fy (2-33) 
A 



MM' 



M T M' 



f d 3 q 1 

x / TTTTT 2 — i ^2 ^ ~ h', J MTMTI70757 • qT\ex.p(iq ■ n)|0) 
7 (2tt)^ m 2 + g 2 

x(p - h; J'M'T'M' T | 7 o757 • <fr_ A exp(-ig • r 2 )|6) 

= ^(-i) j+1 a^TI)(-(|^) 2 ) 

>< f 0^{Erf(*) + E-^(*)} 

X {E ai '- 7i? p J h'^) +^2 a JL' F ph'(Q)}- 
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We shall include the pion form factor to take into account the finite nucleon size 
effect by replacing the pion propagator as^3 

1 1 /A 2 -ml\2 

f + ml ~^ f +ml V ? + A 2 ) ' { ' 

where A represents the momentum cutoff parameter. 

We would like to show here the integrand of the momentum integral of Op-Oh and 
2p-2h matrix elements to see the connection of the pion quantum number in finite 
system and pion momentum in infinite matter. Shown in Fig. [2] is the momentum- 

2 

dependent part of the pion-exchange interaction, m i +q 2 , together with the form 
factor, and the squared particle-hole transition density, q 2 F^ h (q), as functions of 
momentum, q. Here, q 2 is the factor coming from the phase space, and we multiply 
this factor by F^ h (q) shown in Eq. (2.35). 

= {Y a ^ F ^*(l) + Y, a rL^t{<l)} (2-35) 

L L 

x { Y, a L'jFph (?) + E a JL'^ph (9) } ■ 

V L> 

The 0~ component has the largest transition density and peaks at smaller momenta, 
while the large J w components have gradually smaller transition densities and peak 
at larger momenta. To take into account the full strength of the pion-exchange 
interaction in the intermediate-range region, we have to include these higher spin 
pionic states. The CPPRMF model takes into account the pionic strength indicated 
here by J w = 0~. 




q (fm"') 



Fig. 2. Momentum-dependent part of pion exchange interaction, m 2 +q 2 1 together with the form 

factor, and the squared particle- hole transition density, q 2 Fp h (q), as functions of q. The squared 
particle-hole transition densities are shown for various J n by solid curves and the pion exchange 
interaction is shown by a dashed curve for A = 1000 MeV and by a dotted curve for A — 800 
MeV. The parameters used for this presentation are a meson mass, m CT = 850 MeV and the 
w-nucleon coupling constant, = 7.07. 
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§3. Numerical results and discussions 

We apply the RCMF model to the 4 He nucleus as a pilot calculation. As the first 
step, we prepare Op-Oh RMF ground state, (Os^) 4 , by solving the RMF equation 
obtained from the a-u Hamiltonian, H ajUJ . We adjust the a meson mass, m a , and 
the w-nucleon coupling constant, g w , so as to reproduce the root-mean-square (rms) 
matter radius of 4 He, R Tms = 1.488 fm. We use two ranges for Gaussian function with 
bo = 0.7 and 1.4 fm, where the latter range is the dominant range. We keep these 
two ranges during the self-consistent calculation. For the a and oj meson fields, 
we expand them in the same way as the hole state using the two-range Gaussian 
functions, 

/ 5 / 2 / r 2 \ 

u=l,2 V °<» ^ 0V 



where eft denotes a and u meson fields. 

3.1. One-range Gaussian for particle states 

In Fig. [3l we show the total energy of 4 He and various energy components as 
a function of the Gaussian range of particle state, b. We take 2p-2h states up to 
-'max = 6~ of the pionic quantum number. The total energy minimum is achieved at 
around b = 0.8 ~ 0.7 fm, which corresponds to about half the Gaussian range, 1.4 
fm, which is the dominant range of the hole state. The particle state has a spatially 
compact distribution reflecting high-momentum components in the wave function 
due to the pseudo-scalar nature of the pion-exchange interaction, (<7j • q)(ffj -q). This 
important feature has been pointed out in CPPHF 10 ' and CPPRMF^ 1 methods in 
the case of the spherical pion field ansatz, J 77 = 0~ . It is indispensable to have this 
wide variational space for 2p-2h states to take into account properly the pseudo- 
scalar nature of the pion-exchange interaction. Reflecting the pseudo-scalar nature, 
there is a strong correlation between kinetic energy and pion energy. The central 
energy, which is the summation of the scalar and vector potentials, is almost constant, 
which is not shown in Fig. [3l The energy contribution from the nonlinear term of 
a and u> mesons is also almost constant. We show the results of two cases for cutoff 
momentum of the pion form factor. The contribution of the pionic energy is large 
for large cutoff momentum, A. 

In Fig. [U we show the total energy and its components for 4 He as functions of 
the pionic quantum number, J 71 ", for A = 1000 MeV in the left-hand panel and A 
=800 MeV in the right-hand panel. We see that the total energy and its components 
tend to converge as it approaches J w = 6 _ , which corresponds to the single-particle 
orbital angular momentum, L = 7. These high-multipole components are required 
to take into account the full strength of the intermediate-range part of the relative 
distance of two-nucleon pairs due to the pion-exchange interaction for 4 He. The 
range parameter, which makes total energy minimum, changes from 0.83 fm with 
J* = 0" to 0.75 fm with 6". 
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0.4 0.6 0.8 1.0 1.2 0.4 0.6 0.8 1.0 1.2 



Range of the Gaussian b (fm) Range of the Gaussian b (fm) 

Fig. 3. Total energy and various energy components as functions of the Gaussian range of particle 
state, 6, for 4 He with JJ, ax = 6~. We take two cases for the cutoff momentum A — 1000 MeV 
(left figure) and 800 MeV (right figure). The other free parameters, a meson mass and w-nucleon 
coupling constant, are 850 MeV and 7.07, respectively 




Pionic quantum number J* Pionic quantum number J* 



Fig. 4. Total energy and various energy components as functions of pionic quantum number, J", 
for 4 He. The cutoff momenta, A, are taken to be 1000 MeV (left figure) and 800 MeV (right 
figure). The other free parameters, a meson mass and cj-nucleon coupling constant, are 850 
MeV and 7.07, respectively 



3.2. Multirange Gaussian expansion for particle states 

We show in Fig. 0the total energy and the pion and central energies as functions 
of the number of Gaussian functions. We search the energy minimum by changing the 
Gaussian range, 6j, for each step where we add a new Gaussian function. We increase 
progressively a new Gaussian function until the energy convergence is realized. We 
obtain energy convergence by using a relatively small number (~ 6) of Gaussian 
function. This is an advantage of the pionic optimized method. Because the range 
parameter, bi, is a variational parameter, we can minimize the number of range 
parameters in the above method. Almost all range parameters, bi, are found to be 
less than 1.0 fm. In actual calculations, we take 8 ranges, whose values are bi = 
0.4, 0.6, 0.8, 0.9, 1.0, 1.1, 1.4, 2.0 fm. We include mostly small-size Gaussians and 
one large size, b = 2.0, to take care of noncompact components. The pion energy 
increases by around 20% compared with the case of the one-range Gaussian. In this 
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Fig. 5. Convergence of the total energy with the number of Gaussians for 4 He. The pionic quantum 
number is taken up to J^ax = 1 + - The cutoff momentum, A, is taken to be 1000 MeV. Two free 
parameters, a meson mass and cj-nucleon coupling constant, are 850 MeV and 7.07, respectively. 



figure, we take up to J n = 1 + , where channels of particle states, (pi/2) 2 , (P3/2) 2 i an d 
(lsi/2)(^3/2)> are taken into account, where the hole state is (Os^) 2 - Hereafter, we 
do not write explicitly the hole state. 

3.3. Nuclear size and energy 

We show in Fig. [6] the rms matter radius of 4 He as a function of the a meson 
mass in three cases with the cutoff momenta, A = 700, 800, and 1000 MeV. We 
adjust the w-nucleon coupling constant to reproduce the empirical binding energy of 
4 He without the Coulomb interaction. As the cutoff momentum, A, becomes larger, 
a heavier a meson mass is required to obtain the rms matter radius. When the a 
meson mass becomes larger, the attractive force decreases, while the pion-exchange 
interaction contributes more to the binding energy. As a consequence, the equation 
of state becomes softer, as the cutoff momentum becomes larger. In the case of A = 
1000 MeV, the rms matter radius gradually becomes flat as the a meson mass. Even 
if we take a significantly large a mass, it is not possible to find m a to reproduce the 
binding energy and matter radius simultaneously. In cases of A = 700 and 800 MeV, 
the a meson mass under the realization of the saturation property is significantly 
large compared with that (500 ~ 600 MeV) of the Walecka modef^'^ and that 
(850 MeV) obtained under the spherical pion field ansatz, J 71 " = 0~ , in the CPPRMF 
method.^ In the case where the matter radius and total energy are reproduced, we 
find that the bound state cannot be realized only by scalar and vector potentials. The 
bound state is realized by the attraction almost entirely due to the pion-exchange 
interaction. 

In FigJTl we show various energy components of 4 He as functions of the a meson 
mass. As the a meson mass becomes heavier, its field strength decreases and the 
attraction due to the a meson is reduced. Instead, the energy contribution due to the 
pion-exchange interaction increases. The reduction of the strength of the a meson is 
accompanied by a decrease in the strength of the uj meson field. As a consequence, the 
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Fig. 6. Root-mean-square matter radius of 4 He as a function of the a meson mass for three cases 
of the cutoff momentum, A. The result of A — 1000 MeV is shown by the solid circle, that of 
A = 800 MeV by the solid square, and that of A = 700 MeV by the open square. The pionic 
quantum number is taken up to J^ax = 6 _ . The number of Gaussians is 8, where bi = 0.4, 0.6, 
0.8, 0.9, 1.0, 1.1, 1.4, 2.0 fm. 
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Fig. 7. Various energy components for 4 He as functions of the a meson mass, m CT , with the cutoff 
momentum of A = 800 MeV. The u;-nucleon coupling constant, g u , is adjusted to reproduce the 
empirical binding energy of 4 He without the Coulomb interaction. The pionic quantum number 
is taken up to J„ lax = 6~. The number of Gaussians is 8, where 6j = 0.4, 0.6, 0.8, 0.9, 1.0, 1.1, 
1.4, 2.0 fm. 



strength of the spin-orbit force becomes smaller as the a meson mass increases. This 
fact indicates the importance of pions in the nucleus on the formation of the shell 
structure (spin-orbit splitting effect Around 70% of the total attractive 
potential comes from the pion-exchange interaction at the a meson mass, m a = 
1250 MeV in the case of A = 800 MeV. This amount of energy contribution due to 
the pion-exchange interaction is in good agreement with the result of the variational 
Monte-Carlo (VMC) method for light nuclei (A < 10) obtained by the Argonne- 
Illinois group. 16 ' They point out that the energy contribution of the pion-exchange 
interaction to the net two-body attraction is 70~80%. 
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3.4. Tensor and spin-spin central parts of pion- exchange interaction 

We new discuss the contributions of the tensor part and spin-spin central part of 
the pion-exchange interaction. The pion-exchange interaction is separated into the 
spin-spin part and tensor part as 

ai ■ qa 2 ■ qn ■ t 2 l/_ _ q 2 3a x • qa 2 ■ q- q*o\ ■ o 2 \ ,„ ^ 
=2- o = ol 0- ! ' <72 ^2~1 2 H =2~~i 2 Jn-T 2 . {3-2) 

The first term is the spin-spin central part and the second term is the tensor part. 

We show the energy contributions from the tensor part (upper half of left-hand 
panel) and spin-spin part (lower half of left-hand panel) of the pion-exchange in- 
teraction for various 2p-2h states in Fig. [HJ In the right-hand panel, we show 
the energy contributions of the tensor part and spin-spin central part of the pion- 
exchange interaction as functions of J 7r . The most important 2p-2h configuration is 
the (pi/ 2 ) 2 state. This channel corresponds to the case of the spherical pion field 
ansatz, J v = 0~ , which is only the configuration taken into account in the CPPRMF 
modeP 1 '^ 1 ® and also some other works^®'^ in the nonrelativistic version. The 
tensor part behaves similarly to that determined by Myo et alP^f On the other 
hand, we see a large contribution of the (p3/ 2 ) 2 state due to the spin-spin central 
interaction. This spin-spin central contribution being large is not seen in other calcu- 
lations PJ In our results, the tensor and spin-spin central parts give almost the same 
amount of energy contribution for A = 800 MeV as shown in the right-hand panel 
of Fig. [51 About 40% of the total attractive potential comes from the tensor part, 
and about 30% comes from the spin-spin central part. We compare these energy 
contributions with the result obtained by the Argonne- Illinois group P^ 1 The energy 
contribution of the tensor part to the net two-body attractive force is approximately 
50% and that of the spin-spin central part is just approximately 20%. From this 



> 




Channel Pionic quantum number J' 



Fig. 8. Energy contributions of the tensor part and spin-spin central part of the pion-exchange 
interaction for various 2p-2h states. In the left-hand panel shown in the upper half is the tensor 
part and in the lower half is the spin-spin central part obtained with A = 800 MeV. In the right- 
hand panel is shown the energy contributions of the tensor part and spin-spin part of the pion- 
exchange interaction as functions of J*'. The numbers shown in the horizontal axis of the left- 
hand panel represent 2p-2h states; l:(pi/ 2 ) 2 , 2:(ls 1/2 )(d 3/2 ), 3:(d 3/2 ) 2 , 4:(p 3/2 ) 2 , 5:(p 3 / 2 )(/ 5 / 2 ), 
6:(/ 5/2 ) 2 , 7:(d 5/2 ) 2 , 8:(d 5/2 )(<? 7/2 ), 9:( ff7/2 ) 2 , 10:(/ 7/2 ) 2 , ll:(/ r/2 )(ft 9/2 ), 12:(h 9/2 ) 2 , 13:(<? 9/2 ) 2 , 
14:(ff 9 / 2 )(iii/ 2 ), 15:(in/ 2 ) 2 , 16:{h 11/2 ) 2 , 17:(/lu/ 2 )(ii 3 / 2 ), 18:(ji 3/2 ) 2 . The number of Gaussians 
is 8, where b t = 0.4, 0.6, 0.8, 0.9, 1.0, 1.1, 1.4, 2.0 fm. 



16 



Y. Ogawa, H. Toki, S. Tamenaga and A. Haga 



comparison, we can conclude that the energy gain due to the tensor part is under- 
estimated by around 10%; on the other hand, that due to the spin-spin central part 
is overestimated by around 10% in our calculation. 

The reason why we have overestimated the energy due to the spin-spin central 
part is that we do not consider at all the effect of the short-range repulsion, which 
was pointed out by JastrowpSJ which cuts down the wave function in the region 
less than 0.5 fm in the relative distance of two nucleons. The effect of the short- 
range repulsion is serious for the spin-spin central part, because the spin-spin central 
part works when the change in the relative orbital angular momentum between two 
nucleons is 0, AL = 0. Since there is no centrifugal barrier, its amplitude becomes 
large, as two nucleons approach each other. On the other hand, the tensor part 
works when the change in the relative orbital angular momentum is 2, AL = 2. The 
relative wave function of a two-nucleon pair for the tensor case is pushed outside by 
the centrifugal potential. 

We consider that the problem is related to the overestimation of the spin-spin 
central part. We have not treated at all the effect of the short-range repulsion. 
If we perform the correct procedure for taking the many-body correlation in the 
short-range part, we may avoid the overestimation of the spin-spin central part. 
As a suitable procedure, the unitary operator correlation method (UCOM) is well 
knowrP^ 1 to describe the short-range correlation in the framework based on the 
single-particle picture. Once we can suppress the contribution of the spin-spin central 
part, we can increase the contribution of the tensor part by taking a larger cutoff 
momentum, A. We discuss UCOM in the following section. 

3.5. Inclusion of the effect of the short-range repulsion 

We follow the method of UCOM for the spin-spin central part of the pion- 
exchange interaction, which is the first term in Eq. (3.2) as 

V(q) = ^j-^. (3-3) 
q l + 7774 

In UCOM, a unitary correlation operator is introduced as a shift operator, C, where 
■ip = C(fi represents correlated states and C^HCtfi = E(j> represents the equation of 
the uncorrelated states, 4>. This unitary correlation operator moves two particles 
away from each other whenever they are within a short distance due to the repul- 
sive coreP3 We apply this unitary correlation operator for only the central part 
of the pion-exchange interaction within the approximation of taking terms up to 
two-particle correlated operators. The interaction with form factor in coordinate 
representation is 

1 f /A 2 — m 2 \ 2 

v<r) = ToTvT / Whrr5 ^(^f f) (3-4) 



(2tt) 3 J y A 2 + q 2 

1/ 2 exp(-m 7r r) A\ - m 2 ^ expj-Axr) A\ - m 2 A2 exp(-A 2 r) \ 
^{ m * r A 2 -A 2 A r + A 2 - A 2 ^ 2 r J' 
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Fig. 9. Root-mean-square matter radius of 4 He for several cases of cutoff momentum, A as function 
of the a meson mass in the left panel. The solid circle represents the result for A — 1000 MeV, 
the solid square for A — 800 MeV, and the open square for A = 700 MeV. The right-hand panel 
shows energy components as functions of the a meson mass for A = 1000 MeV. The pionic 
quantum number is taken up to J max = 6~ . The number of Gaussians is 8, where bi = 0.4, 0.6, 
0.8, 0.9, 1.0, 1.1, 1.4, 2.0 fm. 



where A\ = A + e, Ai = A — e, and e is taken to be 10 MeV to express differential 
operation with respect to A^S The correlated operator is given as 

V(r) = C^V{r)C = V(R+(r)), (3-5) 

where R+(r) = r + a(^) v exp(— exp(^)). The parameters, a = 0.94 fm, (3 = 1 
fm, and r] = 0.37 are used in Ref. 27). The correlated operator in momentum 
representation, V(q), is obtained by Fourier transformation of V(r). 

We show the mean square radius of 4 He in Fig. To reproduce both the 
binding energy and the rms matter radius simultaneously, we need a significantly 
large a meson mass as compared with the case of the phenomenological parameter 
set, for example, the Walecka modelP^^ As the a meson mass becomes heavier, 
the attraction due to the a field decreases and the energy contribution from the 
pion-exchange interaction increases. The a meson mass should be m a = 1300, 1000, 
and 940 MeV in cases with A = 1000, 800, and 700 MeV, respectively, to reproduce 
both the rms matter radius and binding energy. In the right-hand panel of Fig. [9j 
we show the energy components for 4 He in the case of A = 1000 MeV. Around 70% 
of the total attractive potential comes from the pion-exchange interaction when the 
binding energy and rms matter radius are reproduced simultaneously. This amount 
of energy contribution due to the pion-exchange interaction is in good agreement with 
the result of the VMC method for the light nuclei obtained by the Argonne- Illinois 
group. ™ The energy contribution due to the pion-exchange interaction decreases as 
the cut-off momentum becomes smaller. 

We calculate the energy contributions from both the tensor part and spin-spin 
central part of the pion-exchange interaction given in Eq. (3.2) and show the results 
in Fig. [TOl When we take into account the effect of the short-range repulsion by using 
the UCOM prescription for the spin-spin central part, the ratio of the energy from 
the spin-spin central part to the one from the tensor part is around 1/3. Around 



18 



Y. Ogawa, H. Toki, S. Tamenaga and A. Haga 



I I 



12 13 14 15 16 17 



-14 

12 

_ -10 
> 

g> -6 
W - 4 

-2 




Spin-spin central 



with UCOM 
without UCOM 



12 13 14 15 16 17 18 



Fig. 10. Energy contributions of the tensor part and spin-spin central part of the pion-exchange 
interaction for various 2p-2h states with A — 1000 MeV. The left-hand panel shows the energy 
from the tensor part. In the right-hand figure is shown the energy from the spin-spin central 
part. The thin bars in the right-hand panel represent the energy without the UCOM operation. 
The binding energy and rms matter radius are reproduced simultaneously by adjusting the a 
meson mass and aj-nucleon coupling constant. Each number shown in the horizontal axis rep- 
resents the 2p-2h channels, where l:(pi/ 2 ) 2 , 2:(ls 1/2 )(d 3/2 ), 3:(d 3/2 ) 2 , 4:(p 3/2 ) 2 , 5:(p 3 / 2 )(/ 5/2 ), 
6:(/ 5/2 ) 2 , 7:(d 5/2 ) 2 , 8:(d 5/2 )(g 7/2 ), 9:( 57/2 ) 2 , 10:(/ 7/2 ) 2 , ll:(f 7/2 )(h a/2 ), l2:(h 9/2 ) 2 , 13:(<? 9/2 ) 2 , 
14: (59/2) (in/2), 15:(in/ 2 ) 2 , 16:{h 11/2 ) 2 , 17:(/ln/ 2 )(il3/2), 18:(ji 3/2 ) 2 . The number of Gaussians 
is 8, where bi = 0.4, 0.6, 0.8, 0.9, 1.0, 1.1, 1.4, 2.0 fm. 



50% of the total attractive potential comes from the tensor part, and about 20% 
comes from the spin-spin central part. We compare these energy components with 
the result obtained by the Argonne-Illinois group. 16 ' We obtain good agreement 
with their result. In particular, we can reduce the energy of the spin-spin central 
part of the (P3/2) 2 state. 




01234 01234 
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Fig. 11. The vector density distribution is shown as a function of the radial distance, r, in the left 
panel for A — 1000 MeV and in the right panel for A — 800 MeV. The dash-dotted curve repre- 
sents the density of the Op-Oh RMF ground state, while the dashed curve represents the density 
distribution of the two-particle state and the dotted curve represents the density distribution 
of the two-hole state. The solid curve represents the density distribution of total ground state 
wave function, Op-Oh + 2p-2h. For both cases, the parameters are adjusted to reproduce the 
binding energy and rms matter radius simultaneously. The pionic quantum number is taken up 
to Jm ax = 6". The number of Gaussians is 8, where h = 0.4, 0.6, 0.8, 0.9, 1.0, 1.1, 1.4, 2.0 fm. 
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3.6. Density distribution 

In Fig. [TT1 we show the vector density distribution for 4 He for A = 1000 MeV 
in the left panel and for 800 MeV in the right panel. In both cases, the binding 
energy and rms matter radius are reproduced simultaneously. The dashed and dotted 
curves represent the densities of two-particle component and two-hole component of 
2p-2h states, respectively. The integrated value of the density of the two-particle 
component is twice the probability of the 2p-2h states, X^ a i a «- The amount of 
this probability is around 14% for A = 1000 MeV. The amount of twice of this 
probability of the density distribution moves from the Op-Oh RMF ground state to 
the two-particle excited state through the pionic correlation due to the pion-exchange 
interaction. This 2p-2h contribution changes the density distribution of 4 He. The 
central part of the density distribution is reduced and the outside region at r > 2 fm 
of the density distribution is also reduced, and consequently, the distribution becomes 
spatially compact. A peak of the fluctuation density distribution of the particle-state 
arises at around 0.8 fm, and this fact indicates the important property of the pion- 
exchange interaction. The density distribution of two-particle states is spatially 
compact reflecting the pseudo-scalar nature of the pion-exchange interaction. 

We show in Fig. [12] the effective masses for the nucleon and u meson as a 
function of the nuclear radius for A = 1000 and 800 MeV. The effective masses are 
given for the nucleon and oj meson as 

M* = M + g a a, (3-6) 

The pionic energy contribution becomes important in cases of large cutoff momen- 
tum. For large A, the a meson mass should be larger and the contribution of the a 
meson becomes small. Hence, the nucleon and uj meson masses are larger for larger 
A. 
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Fig. 12. Effective masses of the nucleon and a>-meson for A = 1000 and 800 MeV. For both cases, the 
parameters are adjusted to reproduce the binding energy and rms matter radius simultaneously. 
The pionic quantum number is taken up to J^ax = 6 _ . The number of Gaussians is 8, where 
h = 0.4, 0.6, 0.8, 0.9, 1.0, 1.1, 1.4, 2.0 fm. 
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§4. Summary 

We have formulated the RCMF model for finite nuclei based on the relativistic 
single-particle framework. This RCMF model is a natural extension of the CP- 
PRMF model to include higher pionic quantum numbers of particle-hole states, 
J w = 1 + , 2~ , 3 + , .... , to take into account the full strength of the pion-exchange in- 
teraction. The CPPRMF model corresponds to the J n = 0~ case in the 2p-2h states 
in the RCMF model, which reflects the surface contribution of the pion-exchange in- 
teraction. The inclusion of higher multipoles comes from the necessity to include the 
volume contribution of the pion-exchange interaction. We have applied the RCMF 
model to the 4 He nucleus as a pilot calculation and studied the energy convergence 
and role of pion in the nuclear structure. The RCMF model has three free parame- 
ters in the sigma model Lagrangian. They are the a meson mass, m a , the w-nucleon 
coupling constant, g u , and the cutoff momentum, A. These parameters are fixed so 
as to reproduce the binding energy and radius of He and the contribution of the 
tensor interaction of the pion-exchange interaction. We have obtained the energy 
convergence when we take into account the pionic quantum number up to J max = 
6-. 

The Gaussian ranges are taken as the energy variational parameters. The spa- 
tially compact distribution of two-particle states is obtained explicitly at the point 
where the minimum of the total energy is realized. This fact reflects the involvement 
of the high-momentum component in the wave function. This framework can de- 
scribe the pseudo-scalar nature of the pion-exchange interaction in this variational 
method. When we increase the number of Gaussians to expand particle states, the 
pion energy gain increases by around 20% from the case of one-range Gaussian. In 
this framework, the number of Gaussian expansions does not need to be high, be- 
cause the Gaussian range, bi, is taken as the variational parameter. We can take 
into account the full strength of the pionic correlation in the intermediate range by 
using this framework. 

To reproduce both binding energy and rms matter radius simultaneously, we 
need a significantly large a meson mass as compared with that of the phenomeno- 
logical parameter set, for example, the Walecka modelP^'EB As the a meson mass 
becomes heavier, the attraction due to the a field decreases and the energy contri- 
bution from the pion-exchange interaction increases. The a meson mass is taken to 
be m a = 1300, 1000, and 940 MeV in the case with A = 1000, 800, and 700 MeV, 
respectively. 

Around 70% of the total attractive potential comes from the pion-exchange 
interaction for A = 1000 MeV, and at this moment, the binding energy and rms 
matter radius are reproduced simultaneously. This amount of energy contribution 
due to the pion-exchange interaction is in good agreement with the result of the 
VMC method for the light nuclei obtained by the Argonne-Illinois groupP^ 1 We have 
calculated various energy contributions from both spin-spin central part and tensor 
part of the pion-exchange interaction given in Eq. (3.2). Both parts give almost the 
same amount of energy contribution when we do not take care of the short-range 
repulsion. When we take into account the effect of the short-range repulsion by 
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using the UCOM prescription for the spin-spin central part, the ratio of the energy 
contribution from the spin-spin central part and tensor part is around 1/3. Around 
50% of the total attractive potential comes from the tensor part, and about 20% 
comes from the spin-spin central part. We compare these energy constituents with 
the result obtained by Argonne-Illinois group™ We obtain good agreement with 
their result. 

As for the energy contribution due to the tensor part, the most important two- 
particle state in the 2p-2h channels is the (j>i/2) 2 - The next ones are the (P3/2)(/5/2)> 
(l s i/2)(^3/2)> an d {d 3 / 2 ) 2 in order of importance. As the energy contribution due to 
the spin-spin central part, the (P3/2) 2 state leads to the most important contribution. 
The next ones are the (pi/2) 2 , (^5/2) 2 ! (<^3/2) 2 m order of importance. 

In this work, the lower radial wave function is connected with the upper radial 
wave function with the plane wave relation through ^jj, to minimize the number of 
variational parameters. We would like to remove this constraint for small components 
as variational parameters. Furthermore, we have to include the exchange (Fock)- 
term in this framework, since we describe the problem of the many-body system. 
These points will be improved in forthcoming works. We shall proceed to calculate 
heavier nuclei to ensure the validity of our statement in the previous work, where the 
pion plays a role in the formation of the jj-magic shell structure.® We have shown 
the effective masses of the nucleon and uj meson inside the nucleus. We would like 
to consider the relation between pion and <7-field for the purpose of the discussion 
of the chiral symmetry in finite nucleus in future work. We have applied UCOM 
prescription to treat properly the short-range repulsion for the spin-spin central part 
of the pion-exchange interaction. We have to develop a dynamical framework to 
determine the amount of short-range correlation by explicitly including short-range 
repulsive interaction. 
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